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We derive exact expressions for the mean value of Meyer- Wallach entanglement Q for localized 
random vectors drawn from various ensembles corresponding to different physical situations. For 
vectors localized on a randomly chosen subset of the basis, (Q) tends for large system sizes to a 
constant which depends on the participation ratio, whereas for vectors localized on adjacent basis 
states it goes to zero as a constant over the number of qubits. Applications to many-body systems 
and Anderson localization are discussed. 
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Random quantum states have recently attracted a lot 
of interest due to their relevance to the field of quantum 
information. Since they are useful in various quantum 
protocols [l[ , efficient generation of random and pseudo- 
random vectors |2j and computation of their entangle- 
ment properties |3[ have been widely discussed. 

Random states are not necessarily uniformly spread 
over the whole Hilbert space. It is therefore natural to 
study entanglement properties of random states which 
are restricted to a certain subspace of Hilbert space, or 
whose weight is mainly concentrated on such a subspace. 
Such states can appear naturally as part of a quantum 
algorithm, or can be imposed by the physical implemen- 
tation of qubits, through e. g. the presence of symmetries. 

In addition, random states built from Random Matrix 
Theory (RMT) have been shown to describe many prop- 
erties of complex quantum states of physical systems, es- 
pecially in a regime of quantum chaos. Yet in many cases 
physical systems display wavefunctions which are local- 
ized preferentially on part of the Hilbert space. This 
happens for example if there is a symmetry, or when 
the presence of an interaction dclocalizes independent- 
particle states inside an energy band given by the Fermi 
Golden Rule. A different case concerns Anderson local- 
ization of electrons, a much studied phenomenon where 
wavefunctions of electrons in a random potential are ex- 
ponentially localized. Assessing the entanglement prop- 
erties of such states not only enables to relate the en- 
tanglement to other physical properties, but also has a 
direct relationship with the algorithmic complexity of the 
simulation of such states. Indeed, it has been shown Q 
that weakly entangled states can be efficiently simulated 
on classical computers. 

For a vector in a iV-dimensional Hilbert space, lo- 
calization can be quantified through the Inverse Partici- 
pation Ratio (IPR) £ = Y,i l*i| 2 /Ei l^l 4 wher e are 
the components of This measure gives £ = 1 for a 
basis vector, and £ = M for a vector uniformly spread on 
M basis vectors. 

To investigate entanglement properties of localized 
vectors, we choose the measure of entanglement pro- 
posed in @. Meyer- Wallach entanglement (MWE) Q 
can be seen as an average measure of the bipartite en- 
tanglement (measured by the purity) of one qubit with all 



others. The quantity Q has been widely used as a mea- 
sure of the entangling power of quantum maps or 
to measure entanglement generation in pseudo-random 
operators @. For a pure A^-dimensional state \P coded 

on n qubits [N = 2"), Q = 2 (l - ± £™=o i? r ), where 

R r = tvpl is the purity of the r-th qubit (p r is the par- 
tial trace of the density matrix over all qubits but qubit 
r). It can be rewritten as Q = — E™=o G{u r , v r ), where 
G(u,v) = (u\u)(v\v) — \ (u\v}\ 2 is the Gram determinant 
of u and v, and u r (resp. v r ) is the vector of length N/2 
whose components are the VE^ such that i has no (resp. 
has a) term 2 r in its binary decomposition. Vectors u r 
and v r are therefore a partition of vector 'F in two sub- 
vectors according to the value of the r-th bit of the index. 

Analytical computations will be made on ensembles 
of random vectors. In this case, individual quantum 
states in a given basis have components whose ampli- 
tudes, phases and positions in the basis are drawn from 
a distribution according to some probability law. Quan- 
tities such as IPR or entanglement measure are then av- 
eraged over all realizations of the vector. A simple exam- 
ple of a random vector localized on M basis states can be 
constructed by taking M components with equal ampli- 
tudes and uniformly distributed random phases, and set- 
ting all the others to zero. A more refined example con- 
sists in using, as nonzero components, column vectors of 
M xM random unitary matrices drawn from the Circular 
Unitary Ensemble of random matrices (CUE vectors). 

In the first part of this paper we study entanglement 
properties of random quantum states which are localized, 
or mainly localized, in some subset of the basis vectors. 
We show that very different behaviors can be obtained 
depending on the precise type of localization discussed. 
The first case we consider (section [TJ) consists in random 
states whose non zero components in a given basis are 
randomly distributed among the basis vectors. Moreover, 
these nonzero components are chosen to have random val- 
ues. Averages over random realizations therefore imply 
that we average both over position of the nonzero com- 
ponents among the basis vectors and over the random 
values of these nonzero components. We show that the 
mean entanglement can be expressed as a function of the 
number of nonzero components of the vector. We then 
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show that this result can be generalized. Indeed for any 
vector with random values distributed according to some 
probability distribution, the mean entanglement can in 
fact be expressed as a function of the mean IPR. Notably, 
this function tends to a constant close to 1 for large sys- 
tem size. While the vectors in section U are localized on 
computational basis states which are taken at random, in 
section [TT] random vectors are localized on computational 
basis states which are adjacent when the basis vectors are 
ordered according to the number which labels them. In 
this case the mean entanglement can again be expressed 
as a function of the mean IPR, but in contrast this func- 
tion tends to for large system size. Again, the averages 
are performed both on position and values of the compo- 
nents. In the second part of the paper, we compare these 
results to the entanglement of various physical systems 
which display localization ( section HlT|) . 

The question of entanglement properties of localized 
states has already been addressed in other works. The 
concurrence of certain localized states in quantum maps 
has been studied in [!, [f| , but with an emphasis on the 
effect of noise in quantum algorithms. In [7|, a relation 
between the linear entropy and the IPR has been derived 
in the special case where each qubit is an Anderson lo- 
calized state. During the course of this work, a preprint 
appeared which uses different techniques to relate the en- 
tanglement to the IPR 8] in the case of vectors localized 
on non-adjacent basis states, as in section |TJ Interest- 
ingly enough, the formulas obtained in Q are fairly gen- 
eral. They are derived by different techniques and rest 
on different assumptions. In particular, the authors of [|[ 
do not average over random phases. They obtain a for- 
mula where entanglement is expressed as a function of the 
mean IPR calculated in three different bases, a quantity 
that is often delicate to evaluate. Our work uses differ- 
ent techniques and the additional assumption of random 
phases to get a different formula (formula (3)) which in- 
volves only the IPR in one basis, a quantity that can be 
easily evaluated in many cases and is directly related to 
physical quantities such as the localization length. For 
example, it enables us to compute readily the entangle- 
ment for localized CUE vectors (see (4)). However there 
are instances of systems (e.g. spin systems) where these 
different formulas give the same results. 



I. ANALYTICAL RESULTS FOR RANDOMLY 
DISTRIBUTED LOCALIZED VECTORS 

Let us first consider a random state \I> of length N = 2 n 
in the basis {\i} — \io) ® • • • ® \i n -i), < i < 2 n — l,i = 
Y^r=o v2 r } of register states (where all a* are diagonal). 
Suppose the state \t has M nonzero components which 
we denote by ipi, 1 < i < M. Each nonzero component 
is random and additionally corresponds to a randomly 
chosen position among basis vectors. The corresponding 
average will be denoted by (...). We make the assumption 
that these components have uncorrelated random phases, 



and that (|Vvl 2 ) an d (l^pl 2 !^! 2 ) do not depend on p,q. 
We calculate the contribution to MWE of a partition 
(it, v) (we drop indices r). Suppose u has k non-zero 
components Uj, i € I and that v has M — k non-zero 
components vj, j 6 J, with I, J subsets of {1, N/2}. 
We define T = I D J and the bijections a and r such 



that Ui = ipvU) and Vj 



average G{u,v) is given by 



Setting s p = \ip p \ , the 



V^StW^ (1) 



have vanished 



{G(u,v)) = { s p s y 

pea(I) qer(J) 

where the non-diagonal terms in \(i 
by integration over the random phases of the ip p . We 
assumed that (s p s q ) (p ^ q) does not depend onp, q, thus 
(G(u, v)) = [k(M — k) — t](s p s q ), the overlap t being the 
number of elements of T. Since {u\u} + (v\v) = 1, we also 
have (G(u,v)) = k((s p ) - (s»» - [k(k- 1) +t]{s p s q ). We 
then equate both expressions and use our hypothesis that 
(\ip P \ 2 ) and (|V'p| 4 } are independent of p, which implies 
that (s p ) = l/M and {si) = (l/£)/M, to get 



k(M -k)-t 1 
< G ("'»» = M(M-l) 



(2) 



As this result depends only on (k,t), the calculation of 
(Q) comes down to counting the number of positions of 
the non-zero components in vectors u and v yielding the 
same pair (k,t). The combinatorial weight associated to 
a given {k,t) is (^ 2 ) (J) (^ a fc ~* ) . At fixed k, t ranges 
from to min(fc, M 
yields: 

(Q) = 



k). Summing all contributions 



N -2 ( .1. 



(3) 



This result does not depend on M. It can in fact be 
derived by an alternative method with less restrictive as- 
sumptions. Let us sum up all the localization properties 
of iff in the IPR £ alone, whatever the value of M. We de- 



fine the correlators C xx 

Gxy 



(k 



2 )/2, and 



Uj| 2 |t;j| , where the over line denotes the average 
taken over all n partitions (u r 1 v r ) corresponding to the 
n qubits, and over all i,j € {1, ...,N/2} with i ^ j (for 
C xx ) and alii, j e {1, ...,N/2} (for C xy ). Thus C xx quan- 
tifies the internal correlations inside u and v, and C xy the 
cross correlations between u and v. Normalization im- 
poses that (l/0+N(N/2-l){C xx ) + (N 2 /2)(C xy ) = 1, 
and Eq. Q leads to (Q) = N(N - 2){C xy ). The as- 
sumption (C xx ) = (C X y) is then sufficient to get Eq. ([3]). 
This derivation also shows that if the phases are uncor- 
related and formula ([3]) does not apply, then necessarily 

(C xx ) (C X y). 

Our result Eq. ([3]) involves only the mean IPR in one 
basis, and uses the assumptions that on average cross 
correlations are equal to internal correlations for the par- 
titions, whatever the probability distribution of the com- 
ponents, and that random phases are uncorrelated. This 
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is to be compared with the result in [8| where (Q) is 
related to the sum of IPR for three mutually unbiased 
bases. Their result does not use the assumption of un- 
corrected random phases, but requires a stronger hy- 
pothesis on correlations (namely, that vector component 
correlations in average do not depend on the Hamming 
distance between the corresponding vector component in- 
dices). In particular, our formula ^ allows to compute 
(Q) e. g. for a CUE vector localized on M basis vectors; 
in this case £ = (M + l)/2, and we get 



(Q) = 



M — IN —2 
M + 1 N - 1 ' 



(4) 



In [llj , Lubkin derived an expression for the mean MWE 
for non-localized CUE vectors of length N, giving (Q) 
(N — 2)/(N + 1). Consistently, our formula yields the 
same result if we take M — N. For a vector with constant 
amplitudes and random phases on M basis vectors, £ = 
M and 



(Q) 



M -IN -2 
M N -Y 



Formula ((3]) can be easily modified to account for the 
presence of symmetries. For instance, suppose the sys- 
tem presents a symmetry which does not mix basis states 
within two separate subspaces of dimension N/2. It is 
then easy to check that N in ([3|) should be replaced by 
N/2. 



II. ANALYTICAL RESULTS FOR ADJACENT 
LOCALIZED VECTORS 

Up to now we have considered random vectors whose 
components were distributed over a randomly chosen 
subset of basis vectors. However in many physical sit- 
uations vectors are localized preferentially on particular 
subspaces of Hilbert space. An important case consists 
in random vectors localized on M computational basis 
states which are adjacent when the basis vectors are or- 
dered according to the number which labels them. The 
general form of such a vector would be \c),...,\c + M — 1), 
< c < 2 n — 1. Again, averaging over random realiza- 
tions of the coefficients of 'J we get Eq.®. The calcu- 
lation of (Q) therefore reduces to determine k and t for 
all qubits and all possible choices of the basis vectors on 
which "J has nonzero components. For a given r, vectors 
u and v correspond to a partition of the set of the com- 
ponents ^ of \& according to the value of the r-th bit of 
i. For instance for the qubit r = 1, and M — 9, N = 16, 
a typical realization of vectors u and v would be 



u= ( Vi 

V = ( l\)2 "03 



1p4 1p5 1p8 
i>6 i>7 o 



^9 




(6) 



Each vector u and v can be split into 2 n ~ 1 ~ r blocks of 
length 2 r . There are Nn ways of constructing such pairs 
(it, v), by choosing a qubit r and a position c for ipi. The 



numbers k and t depend on three quantities: the label 
r G {0, ...,n — 1} of the qubit whose contribution is 
considered; the position c r G {0, . . . , 2 r — 1} of ipi within 
a block, either in u or in v; the remainder m r of M mod 
2 r+1 . Let r be such that 2 r °- 1 < M < 2 r « . One has 
to distinguish the contributions coming from qubits such 
that < r < ro and qubits such that r > tq. First 
consider < r < ro. Suppose ipi is a component of 
vector u. One can check that I U J has k + t + c r = M 
elements, and I\T has k—t = g r (m r +c r ) elements, where 
g r (x) = 2 r g(x/2 T ) with g(x) = |1 - |1 - x\\, x G [0,3[. 
These two equations lead to k — h(M — c r + g r (m r + c r )) 

Similarly, when ipi is a 
i (M+ c r -g r (m r +c r )) 



and t = \{M — c r — g r (m r + c r )). 
component of vector v, we get k = 



= andi = |(M- 



2 r+1 + g r (m r + Cr))- Altogether this 



(5) (Q) = 



leads to 2 x 2 r different contributions with multiplicity 
2 n_1_r (the number of blocks). If r > ro, t is always zero 
and as the position c r is varied, k runs over {1, M— 1}. 
Summing all contributions together we get 

4(M 
3 _ 



where Xr( x ) 




1)(2" - 2 r °) 



2"+ro 



1 



(7) 



l)/2 r for0< 



x < 2 r . Equation (0 is an exact formula for M < N/2. 
For fixed M and n — > oo, n{Q) converges to a constant 
C which is a function of M and £. For M = 2 r °, ro < n, 



all remainders m. 
to 



(Q) - 



, r < ro are zero, and Eq. |7]) simplifies 
(r + l)M 2 - 2(r Q - 1)M - f 



M{M - I) 



4(M + 1) 
3iV 



1 



1 



(8) 



Numerically, this expression with ro = log 2 M gives a 
very good approximation to Eq. for all M. 

Equation |(7J) is exact for e. g. uniform and CUE vec- 
tors, and we will see in section IIIII that it can be applied 
even when 4" is not strictly zero outside a M-dimensional 
subspace. 



III. APPLICATION TO PHYSICAL SYSTEMS 

We now turn to the application of these results to phys- 
ical systems. Localized vectors randomly distributed over 
the basis states may model eigenstates of a many-body 
Hamiltonian with disorder and interaction. Indeed, the 
latter generically display a delocalization in energy char- 
acterized by RMT statistics of eigenvalues within a cer- 
tain energy range, whereas the distribution of eigenvector 
components is Lorentzian or Gaussian. As an example 
we choose the system governed by the Hamiltonian 



H 



»0i 



2^ ■'•j r! '' J j 

i<j 



(9) 
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This model was introduced in [12| to describe a quan- 
tum computer in presence of static disorder. Here the 
a, are the Pauli matrices for the qubit i. The energy 
spacing between the two states of qubit i is given by 
2Ti randomly and uniformly distributed in the interval 
[Ao — 5/2, Ao + 5/2], The couplings Jy represent a ran- 
dom static interaction between qubits and are uniformly 
distributed in the interval [—J, J]. For increasing interac- 
tion strength J eigenstates are more and more delocalized 
in the basis of register states, and a transition towards 
a regime of quantum chaos takes place, with eigenvalues 
statistics close to the ones of RMT [jjj ■ In parallel, this 
process leads to an increase of bipartite entanglement in 
the system (l3l |. In the following, results will be aver- 
aged over random realizations of the I\ and Jy in ([9]) 
("disorder realizations"), which will be denoted by (...). 

The Hamiltonian ([9]) presents a symmetry which does 
not mix basis states having even and odd number of 
qubits in the state |1). Each symmetric subspace con- 
tains N/2 basis vectors among which for each qubit A/4 
have value [1} and N/4 have value |0). In this case, as 
explained at the end of Section I, N has to be replaced 
by N/2 in |3|). This symmetry has the additional effect 
of making the second term in ([T]) vanish identically for 
all eigenvectors. 

Before applying the results of section [I] to the more 
generic case 5 ~ Ao, we first briefly discuss the specific 
case i < Ao. In this case, the energy spectrum of the 
system is divided into bands corresponding to register 
states with the same number rib of qubits in the |1) state. 
Derealization takes place inside each band separately, 
corresponding to a reduced Hilbert space of dimension 
Nb = (,"J- In this case, all the basis states on which 
the derealization takes place have rib qubits among n in 
the state |1). This implies that the components of the 
wave function are not symmetrically distributed on the 
two vectors u and v of Eq.(TI|. Thus, the correlation as- 
sumption breaks down and Formula ([3]) does not apply. 
However we can derive a specific formula in this case, 
starting back from Eq. ([1]) . The probabilities that a basis 
vector with rib qubits in 1 1) enters into v and u are respec- 
tively r\ = rib/n and 1 — r\. So we expect the norm of u to 
be on average 1 — 77, and the norm of v to be on average r\. 
This implies that for a homogeneously delocalized vector 
one has (Q) — > 477(1 —77) for n — > 00 and 77 constant, since 
the presence of the symmetry makes the second term in 
(TTJ) vanish. Applying this latter formula to the specific 
case rib = 1, we recover the result derived in (tJ. Thus 
(Q) tends to a value between and 1 depending on the 
band, as can be seen numerically in Fig. [TJ 

In the case where 5 m Ao, the bands become mixed 
by the interaction, and derealization takes place inside 
the whole Hilbert space. Formula ^ should apply, once 
modified to take into account the symmetry of the Hamil- 
tonian ([5]) . A straightforward modification of the reason- 
ing leading to Eq. © yields (Q) = (l - (±)). It 

turns out that the presence of this particular symmetry 
allowed the authors of Q to make explicit their formula 



e 




FIG. 1: (Color online) Scaled entanglement with respect to 
the reduced localization length (A^/^)" 1 with Nb = (™)- 
The blue thick solid curves correspond to the second band 
(for n — 11-20), the red thin dashed curves to the third band 
(for n = 11-20), the green thin dotted curves to the fourth 
band (for n = 13, 15, 17), the black thin solid curves to the 
fifth band (for n = 10,11,13,15), the yellow thick dashed 
curve to the sixth band (for n = 12), and the violet thick 
dotted curve to the seventh band (for n = 14). 
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FIG. 2: (Color online) Scaled mean MWE (Q)(N - 2)/N 
of © vs IPR for 5 = A , n = 10 (blue circles) and n = 11 
(green squares). Average is over N/16 central eigenstates and 
100 — 200 disorder realizations. Red solid line is the theory, 
and stars are the data for n = 10 with random shuffling of 
components. Inset: scaled correlator (N/2)(N/2 — l){C X x) 
with same parameters; red line is the result when (C xx ) = 

{C X y). 



in a similar case, yielding the same expression as ours. 

Figure [2] shows the entanglement of eigenvectors of 
Hamiltonian ([9]) compared to this formula. The entangle- 
ment goes to one, but departs from the formula at some 
values of the IPR £. The inset illustrates that this dis- 
crepancy corresponds to a breakdown of the hypothesis 
(Cxx) = (Cxy), because of correlations. These correla- 
tions are probably due to the perturbative regime where 
derealization takes place on a strongly correlated subset 
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FIG. 3: (Color online) Mean MWE vs the number of qubits 
for ID Anderson model with disorder from top to bottom w = 
0.2 (blue), 0.5 (red), 1.0 (green), 1.5 (magenta), 2.0 (cyan), 
and 2.5 (orange). Average is over 10 central eigenstates for 
1000 disorder realizations. Solid lines are the C/n fits of the 
tails. Inset: Value of C = lim n ^oo n(Q) as a function of IPR 
£ (green dots) for the values of w above and w = 0.4, together 
with analytical result from Eq. (0 (red line, top) and from 
Eq. JT0J (blue line, bottom) both for M = 2£ . Stars are the 
C values resulting from a C/n fit of the numerical data for 
CUE vectors of size N with exponential envelope exp(— x/l). 



of states. Figure [2] shows that if these correlations are de- 
stroyed by random permutations of the components, the 
results are in perfect agreement with the theory, even- 
though the distribution of the component amplitudes is 
left unchanged. This confirms that (J3|) can be applied 
if correlations are weak between the vector components, 
whatever their distribution. 

In the case of localization on adjacent basis vectors, 
formula (O can be compared to wavefunctions of elec- 
trons in the regime of Anderson localization. Indeed, 
one dimensional disordered Anderson model is known to 
display localized eigenstates for any strength of disorder. 
This type of localization is a one-body phenomenon, but 
it has been shown that it can be efficiently simulated 
on a n-qubit quantum computer, \1/ describing the par- 
ticle in the position representation [l4j]. The localization 
of the particle takes place on a certain number of adja- 
cent computational basis vectors, and the entanglement 
of the quantum state is related to the entanglement pro- 
duced by the quantum algorithm. The wavefunctions of 
the system are known to have an envelope of the form 
exp(— x/l) where I is the localization length. For N- 
dimensional CUE vectors with such an exponential enve- 
lope, we checked that (Q) is in excellent agreement with 
© with £ = I and M = 2£ (stars in inset of Fig. g]). To 
test the formula on actual wavefunctions of the Ander- 
son model, we consider a one dimensional chain of ver- 
tices with nearest-neighbor coupling and randomly dis- 
tributed on-site disorder, described by the Hamiltonian 
Ha + V. Here Hq is a diagonal operator whose elements 



£i are Gaussian random variables with variance w , and 
V is a tridiagonal matrix with non-zero elements only on 
the first diagonals, equal to the coupling strength, set to 
1. For this system, (...) therefore means averaging over 
the diagonal random values. Figure [3] displays (Q) cal- 
culated numerically for eigenvectors of this system, as a 
function of the number n of qubits for various strengths 
of the disorder w. The expected decrease as C/n is per- 
fectly reproduced for large enough values of n. The inset 
shows the value of the constant C compared to the the- 
ory 0, as a function of £. The deviation from ([7]), in 
particular the saturation for large £, can be understood 
by looking at the structure of eigenvectors in Anderson 
model: when there is no disorder (w = 0) the eigen- 
values are Ek = 2 cos 2iri>k and eigenvectors are plane 
waves with frequency . For weak disorder eigenvectors 
are exponentially localized with localization length £ but 
still oscillate at frequencies distributed as a Lorentzian of 
width 1 /£ around Vk ■ We chose eigenvectors with energy 
Ek ~ {vk « 1/4), yielding rapid oscillations of period 
4 which strongly decrease entanglement. It is easy to 
adapt the analysis leading to Eq. ([7]) for ^ chosen as e. 
g. a vector with ^ j = cos7rj/2, c + 1 < j < c + M, and 
zero elsewhere. For instance for M = 2 r ° , tq < n, we get 
(averaging over c) 



(Q) 



26 4 
~9~ ~ M 



8(3r + l) 4(Af 2 -4)\l 



9M 2 



3M2™ 



(10) 



Asymptotically n(Q) converges to a constant indepen- 
dent of £ = M/2. The inset of Fig. [3] shows that this 
theory captures the behavior of the numerical (Q), al- 
though the saturation constant is different. 
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FIG. 4: (Color online) Mean MWE (left) and IPR (right) 
vs number of qubits for quantum smallworld networks with 
w = 1 and p = 0.001 (blue +), 0.005 (red x), 0.01 (green o), 
0.03 (magenta □), 0.06 (orange *). Logarithm is decimal. 

Let us now add to this system pN links between ran- 
domly chosen vertices. This additional long-range in- 
teraction between few vertices turns the system into a 
quantum smallworld network. Such systems can be ef- 
ficiently simulated on a quantum computer, and display 
a localization-delocalization transition for fixed w when 
p is increased [l5| . Figure Q] shows that this transition 
can be probed through the entanglement of the system. 
Indeed, for small p all eigenstates are exponentially lo- 
calized; (Q) is given by (J7J) and decreases asymptotically 
as 1/n; when p is increased the derealization transition 
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takes place and (Q) is now given by Eq. ([3]): for large n, 
it saturates at 1 — 

In conclusion, we have shown that in localized random 
states the mean MWE can be directly related to the IPR 
£. Entanglement properties are very different if the lo- 
calization is on adjacent basis vectors or not. Compari- 
son with physical systems show that global entanglement 



properties are reproduced, although some discrepancies 
show that they are much more sensitive than e.g. level 
statistics to the details of the system. 
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